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I. INTRODUCTION 

A primary motivation for studying diboson physics is 
that the production of two weak bosons and their in- 
teractions provide tests of the clcctroweak sector of the 
standard model (SM) arising from the vertices involving 
trilinear gauge boson couplings (TGCs) Any devia- 
tion of TGCs from their predicted SM values would be an 



indication for new physics Q and could provide informa- 
tion on a mechanism for electroweak symmetry breaking 
(EWSB). 

The TGCs involving the W boson have been previ- 
ously probed in WW.Wj and WZ production at the 
Tevatron pp Collider [1, 0, H, @ and WW production at 
the CERN e+e~ collider (LEP) 0,1,1 03, at different 
center-of-mass energies and luminosities but no devia- 
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tion from the SM predictions has been observed. The 
LEP experiments benefit from the full reconstruction of 
event kinematics in e + e~ collisions, high signal selection 
efficiencies and small background contamination. At the 
Tcvatron, despite larger backgrounds and limited ability 
to fully reconstruct event kinematics, larger collision en- 
ergies are probed and W Z production can be used to di- 
rectly probe the ZWW coupling. The study of WW and 
WZ production at hadron colliders has focused primarily 
on the purely leptonic final states [H, 0, EH • In this paper 
we present a measurement of the jWW/ZWW couplings 
based on the same dataset used to obtain the recent ev- 
idence for semileptonic decays of WW/WZ boson pairs 
in hadron collisions 

As shown in the tree-level diagrams of Fig. [TJ TGCs 
contribute to WW/WZ production via s-channel dia- 
grams. Production of WW via the s-channel process 
contains both trilinear jWW and ZWW gauge boson 
vertices. On the other hand, WZ production is sensitive 
exclusively to the ZWW vertex. 




FIG. 1: Tree-level Feynman diagrams for the processes of 
WW/WZ production at the Tevatron collider via (a) t- 
channel exchange and (b) and (c) s-channel. 



II. PHENOMENOLOGY 

Unraveling the origins of EWSB and the mass generation 
mechanism are currently the highest priorities in particle 
physics. The SM introduces an effective Higgs potential 
with an upper limit on the Higgs boson mass of ~ 1 TeV 
to prevent tree-level unitarity violation [l3l |. 



In a Higgs-less scenario or for heavier Higgs boson 
masses this unitarity limit on the Higgs boson mass indi- 
cates the mass scale at which the SM must be superseded 
by new physics in order to restore unitarity at TeV ener- 
gies. In this case, the SM is considered to be a low-energy 
approximation of a general theory. Conversely, if a light 
Higgs boson exists, the SM may nevertheless be incom- 
plete and new physics could appear at higher energies. 

The effects of this general theory can be described 
by an effective Lagrangian, £ e ff, describing low-energy 
interactions of the new physics at higher energies in 
a model- independent manner. Expanding in powers of 
(1/Ajvp) Q: 

n>l i NP 

where C^M - m tne gjj(2) L x U{l)y gauge-invariant SM 
Lagrangian, A^vp is the energy scale of the new physics 
and i sums over all operators Oi of the given energy 
dimension (n + 4). The coefficients fi parametrize all 
possible interactions at low energies. Effects of the new 
physics may not be directly observable because the scale 
of the new physics is above the energies currently exper- 
imentally accessible. However, there could be indirect 
consequences with measurable effects; for example, on 
gauge boson interactions. 

For the study of gauge boson interactions, the rele- 
vant terms in Eq. [T] are those that produce vertices with 
three or four gauge bosons. The effective Lagrangian, 
£ c g, that parametrizes the most general Lorentz invari- 
ant VWW vertices (V = Z, 7) involving two W bosons 
can be defined as |15j 

r VWW 

£<=££ — = igY(Wl v W> i V v -lf,W) 

gvww 3 i 1 M" M " 7 

+ iK V WlW v V» v + iJfr-wltWjfV* 
- gYWlW^V" + d»V») (2) 
+ g¥e^P{W;d x W v - dxWlW„)V p 
+ iRyWlWvV^ +i^rWl^W£V' jX 

where e pv \ p is the fully antisymmetric e tensor, W de- 
notes the W boson field, V denotes the photon or Z boson 
field, V„ v = d,y v - d v V^ WV = d„W v - cW M , %u = 
l/2(e fl „xpV Xp ), g 7 ww = —e, and gzww = -ecot 6w, 
where e is the electron electric charge, 9w is the weak 
mixing angle and M\y is the W boson mass. The 14 cou- 
pling parameters of VWW vertices are grouped accord- 
ing to the symmetry properties of their corresponding op- 
erators: C (charge conjugation) and P (parity) conserv- 
ing (gY , Ky, and Ay), C and P violating but CP con- 
serving (g¥), and CP violating (gY , ky, and Ay). In the 
SM all couplings vanish (gY = gY — Ky = Ay = Xy — 0) 
except gY — kv — 1- The value of g\ is fixed by elec- 
tromagnetic gauge invariance (gj = 1) while the value 
of gf may differ from its SM value. Considering the C 
and P conserving couplings only, five couplings remain, 
and their deviations from the SM values are denoted as 
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the anomalous TGCs Agf = (gf — 1), Ak 7 = (k 7 — 1), 
Ak z = (k z — 1), A 7 and Xz- 

If non zero anomalous TGCs are introduced in Eq. [2j 
an unphysical increase in the WW and WZ production 
cross sections will result as the center-of-mass energy, 
y/S, of the partonic constituents approaches Anp- Such 
divergences would violate unitarity, but can be controlled 
by introducing a form factor for which the anomalous 
coupling vanishes as s — > oo: 



Aa(s) 



Aa 



(i + s/A 2 NP r 



(3) 



where n = 2 for 7WW and ZWW couplings, and 
ao is a low-energy approximation of the coupling a(s). 
Thus, the previously described anomalous TGCs scale 
as Aao in Eq. [3] The values of Aao (and ao) are con- 
strained by requiring the 5-matrix unitarity condition 
that bounds the J = 1 partial-wave amplitude of inelas- 
tic vector boson scattering by a constant. These con- 
stants were derived by Baur and Zeppenfeld [l6[ for each 
coupling that contributes to reduced helicity amplitudes 
in WZ, jW, or WW production via s-channel. Calcu- 
lated with M w = 80 GeV, M z = 91.1 GeV and with 
the dipole form factor as given by Eq. [3] the unitarity 
bounds for Ak 7 , Akz, Agf and A TGCs are 



|A«°| < 

IA4I < 
|A<?f°| < 



n n 


1.81 ToV 2 


(„-!)- 




n n 


0.83 TcV 2 


(n-l)« 


-1 A 2 , 

»P , 


n n 


0.84 TeV 2 


- (n-iy 


^NP 



, |AA 7 |< 



0.96 ToV 
=T 

I a \ I n" 0.52 TcV 

|AA Z | S („_l)n-l 



(«-!)' 



J1 JVP 



(4) 

For n = 2 and A^v p = 2 TeV, the unitarity condition sets 
constraints on the TGCs of |A/t°| < 1.81, |AA°| < 0.96, 
| A/4 1 < 0.83, |AA|| < 0.52, and |A 5 f°| < 0.84. The 
scale of new physics, Anp, was chosen such that the uni- 
tarity limits are close to, but no tighter than, the cou- 
pling limits set by data. Clearly, as Anp increases the 
effects on anomalous TGCs decrease and their observa- 
tion requires either more precise measurements or higher 
I. 



III. RELATIONS BETWEEN COUPLINGS 

The interpretation of the effective Lagrangian [Eq. [T] de- 
pends on the specified symmetry and the particle content 
of the underlying low-energy theory. In general, £ e ff can 
be expressed using either the linear or nonlinear realiza- 
tion of the SU(2) l x U(1)y symmetry [17| to prevent uni- 
tarity violation, depending on its particle content. Thus, 
C e g can be rewritten in a form that includes the opera- 
tors that describe interactions involving additional gauge 
bosons, and/or Goldstone bosons, and/or the Higgs field 
and operators of interest for any new physics effects. The 
number of operators can be reduced by considering their 
detectable contribution to the measured coupling. 



Assuming the existence of a light Higgs boson, the low- 
energy spectrum is augmented by the Higgs doublet field 
(j), and SU{2)l and U{X)y gauge fields. Because exper- 
imental evidence is consistent with the existence of an 
SU{2)l x U(1)y gauge symmetry, it is reasonable to re- 
quire £ e ff to be invariant with respect to this symmetry. 
Thus, the second term in Eq. [T] consisting of operators 
up to energy dimension six, is also required to have lo- 
cal SU(2)l x U(1)y gauge symmetry and the underly- 
ing physics is described using a linear realization 18] of 
the SU(2) l x U(\)y symmetry. By considering operators 
that give rise to nonstandard jWW and ZWW couplings 
at the tree level. £ R ff can be parametrized in terms of the 
ai parameters [19(. Those parameters relate to the fi 
parameters of the Lagrangian given in Eq. [T| an d to the 
TGCs in the Lagrangian of Eq. [5] as follows |20l |: 



A4 

Ak 7 = (fw<j> + Jba) o IT" = aw 4> + aBrf, 

Wr - = Ak z + £JFA« 7 = ^ 

Mi 



Agf = fw<t> 



2A a 
N P 



(5) 



A = A 7 = Xz = 3g* 



2Xi 



WWW 



a w 



where g is the SU(2)l gauge coupling constant (g — 
e/smffw), cw — cos 9w, sw = sin 0w, and indices W(j) 
(B(f>) and W refer to operators that describe the inter- 
actions between the W (B) gauge boson field and the 
Higgs field 4>, and the gauge boson field interactions, re- 
spectively. The relations in Eq. [5] give the expected or- 
der of magnitude for TGCs to be 0(M w /A 2 NP ). Thus, 
for Anp ~ 2 TcV, the expected order of magnitude 
for Ak 7 , A<?f , and A is 0(1O -3 ). This gauge-invariant 
parametrization, also used at LEP, gives the following 
relations between the Ak 7 , Agf and A couplings: 

Akz = Agi — Ak 7 • tan 2 9w and A = Xz = A 7 (6) 

Hereafter we will refer to this relationship as the "LEP 
parametrization" [or SU(2)xU(l) respecting scenario] 
with three different parameters: Ak 7 , A and Agf . The 
coupling Akz can be expressed via the relation given 
byEq.El 

A second interpretive scenario, referred to as the equal 
couplings (or ZWW—^WW) scenario [l[, specifies the 
jWW and ZWW couplings to be equal. This is also 
relevant for studying interference effects between the 
photon and iT-exchange diagrams in WW production 
(see Fig. [1]). In this case, electromagnetic gauge in- 
variance forbids any deviation of gj from its SM value 
(Agf = AgJ — ) and the relations between the cou- 
plings become 



An = Akz — Ak 7 and A = Xz = A 7 



(7) 



As already stated, for WW and WZ production the 
anomalous couplings contribute to the total cross section 
via the s-channel diagram. Anomalous couplings enter 
the differential production cross sections through differ- 
ent helicity amplitudes that depend on s. The coupling 
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A primarily affects transversely polarized gauge bosons, 
which is the main contribution to the total cross section. 
Consequently, for a given s, the sensitivity to the cou- 
pling A is higher than to n because A is multiplied by s 
in dominating amplitudes for WW and WZ production. 
Different sensitivity to the n couplings is expected due to 
the choice of scenario: the sensitivity to the k coupling 
in the equal couplings scenario is higher than in the LEP 
parametrization scenario simply because of the different 
relations between Eq. [6] and Eq. [7l 



IV. DO DETECTOR 

The analyzed data were produced in pp collisions at 
y/s = 1.96 TeV by the Tevatron collider at Fermilab 
and collected by the DO detector [U during 2002 - 2006. 
They correspond to 1.07 ± 0.07 fb _1 of integrated lu- 
minosity for each of the two lepton channels (evqq and 
[ivqq). 

The DO detector is a general purpose collider detec- 
tor consisting of a central tracking system, a calorimeter 
system, and an outer muon system. The central track- 
ing system consists of a silicon microstrip tracker and 
a central fiber tracker, both located within a 2 T super- 
conducting solenoidal magnet, with designs optimized for 
tracking and vertexing at pseudorapidities 2^] \t]\ < 3 
and < 2.5, respectively. A liquid-argon and uranium 
calorimeter has a central section covering pseudorapidi- 
ties | r\ | up to k, 1.1, and two end calorimeters that extend 
coverage to [ 77 1 « 4.2, with all three housed in separate 
cryostats [23j . An outer muon system, covering \f]\ < 2, 
consists of a layer of tracking detectors and scintillation 
trigger counters in front of 1.8 T iron toroids, followed 
by two similar layers after the toroids (24[. 

Jets at DO are reconstructed using the Run II cone 
algorithm [25[ with cone radius R — (Ay) 2 + (A0) 2 = 
0.5; where y is the rapidity. Jet energies are corrected 
to the particle level. The jet energy resolution for data, 
defined as a PT /pr, ranges from ~ 15% — 25% for jets 
with p T = 20 GeV to ~ 7% - 12% for jets with p T = 
300 GeV, depending on the rapidity of the jet. 

The DO detector uses a three-level trigger system for 
quickly filtering events from a rate of 1.7 MHz down to 
around 100 Hz that are stored for analysis. Events ana- 
lyzed in the electron channel had to pass a trigger based 
on a single electron or electron+jet(s) requirement, re- 
sulting in an efficiency of 98^3%. The triggers based 
on specific single muon and muon+jet(s) requirements 
are about 70% efficient. Thus, all available triggers were 
used for the muon channel to achieve higher efficiency. 
We select all events that satisfy our kinematic selection 
requirements with no specific trigger requirement. The 
efficiency in this kinematic region is very nearly 100%. 
To estimate and account for possible biases on the shape 
of kinematic distributions, we compare data selected with 
the inclusive triggers to data selected with triggers based 
on a single muon. In the kinematic region of interest, the 



inclusive trigger is estimated to have a shape uncertainty 
of less than 5% and a normalization uncertainty of 2%. 



V. EVENT SELECTION AND CROSS SECTION 
MEASUREMENT 

The analysis presented here builds upon a previous 
publication in which we reported the first evidence of 
WW/WZ production with semileptonic final states at 
a hadron collider [l2T |. Such events have two energetic 
jets from the hadronic decay of either a W or Z bo- 
son as well as an energetic charged lepton and significant 
missing transverse energy (indicating a neutrino) from 
the leptonic decay of the associated W boson. There- 
fore, at the analysis level, we selected events with a re- 
constructed electron or muon with transverse momen- 
tum px > 20 GeV and pseudorapidity \q\ < 1.1 (2.0) 
for electrons (muons), a missing transverse energy of 
$ T > 20 GeV, and at least two jets with p T > 20 GeV 
and 1 77 1 < 2.5. The jet of highest px was required to sat- 
isfy pr > 30 GeV. To reduce background from processes 
that do not contain W— > iu, we required the transverse 
mass [HI from the lepton and $ T to be Mjf > 35 GeV. 
The multijet background, for which a jet is misidenti- 
fied as a lepton, was estimated using independent data 
samples. 

Signal (WW and WZ) and background (W+jets, 
Z+jets, tt and single top quark) processes were mod- 
eled using Monte Carlo (MC) simulation. All MC sam- 
ples were normalized using next-to- leading-order (NLO) 
or next-to-next-to-leading-order predictions for SM cross 
sections, except the dominant background VF+jets, 
which was scaled to match the data as described below. 

In the previously published cross section measurement 
analysis [l2|]> the signal and backgrounds were further 
separated using a multivariate classifier to combine in- 
formation from several kinematic variables. The multi- 
variate classifier chosen was a random forest (RF) classi- 
fier [27], Hll . Thirteen well-modeled kinematic variables 
that demonstrated a difference in probability density be- 
tween signal and at least one of the backgrounds were 
used as inputs to the RF. The effects of systematic uncer- 
tainties on the normalization and on the shape of the RF 
distributions were evaluated for signal and backgrounds. 

The signal cross section was determined from a fit of 
signal and background RF output distributions to the 
data by minimizing a Poisson \ 2 function (i.e., a neg- 
ative log likelihood) with respect to variations of the 
systematic uncertainties [29(, assuming SM jWW and 
ZWW couplings. The fit simultaneously varied the 
WW/WZ and H^+jets contributions, thereby also deter- 
mining the normalization factor for the H^+jets MC sam- 
ple. The measured yields for signal and each background 
are given in Table U and the dijet mass peak extracted 
from data compared to the WW/WZ MC prediction is 
shown in Fig. O The combined fit of both channels to 
the RF output resulted in a measured cross section of 
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20.2 ± 2.5(stat) ± 3.6(syst) ± 1.2(lumi) pb, which is 
consistent with the NLO SM predicted cross section of 
a(WW + WZ) = 16.1 ± 0.9 pb 



TABLE I: Measured number of events for signal and each 
background after the combined fit (with total uncertainties 
determined from the fit) and the number observed in data. 





euqq channel 


[ivqq channel 


Diboson signal 


436 ± 36 


527 ± 43 


VF+jets 


10100 ± 500 


11910 ± 590 


Z+jets 


387 ± 61 


1180 ± 180 


tt + single top 


436 ± 57 


426 ± 54 


Multijet 


1100 ± 200 


328 ± 83 


Total predicted 


12460 ± 550 


14370 ± 620 


Data 


12473 


14392 
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FIG. 2: A comparison of the extracted signal (filled his- 
togram) to background-subtracted data (points), along with 
the ±1 standard deviation (s.d.) systematic uncertainty on 
the background. The residual distance between the data 
points and the extracted signal, divided by the total uncer- 
tainty, is given at the bottom. 



VI. SENSITIVITY TO ANOMALOUS 
COUPLINGS 

For TGCs analysis we use the same selection and set 
limits on anomalous TGCs using a kinematic variable 
that is highly sensitive to the effects of deviations of 
Ak, A, and Agf . Because TGCs introduce terms in the 
Lagrangian that are proportional to the momentum of 
the weak boson, the differential and the total cross sec- 
tions will deviate from the SM prediction in the presence 
of anomalous couplings. This behavior is also expected 
at large production angles of a weak boson. Thus, the 
weak boson transverse momentum spectrum, px, is sen- 
sitive to anomalous couplings and can show a significant 
enhancement at high values oipT- 

The predicted WW and W Z production cross sections 
in the presence of anomalous TGCs are generated with 



the leading order (LO) MC generator of Hagiwara, Zep- 
penfeld, and Woodside (HZW) Q with CTEQ5L [3l| 
parton distribution functions (PDFs). For example, the 
predicted "anomalous" cross sections relative to the SM 
value given by the HZW generator are shown in Fig. [3] as 
a function of anomalous couplings. For this figure we vary 
only the Ak coupling with the constraint between Ak 7 
and Anz as given by Eq. [5] The couplings A and Agf 
are fixed to their SM values (i.e., A = Agf = 0). The 
effects of anomalous couplings on two WW kinematic 
distributions (pt and rapidity of the qq system) for the 
LEP parametrization are shown in Fig. 2J Here again, 
we vary only one coupling at a time (Ak, A or Agf ) ac- 
cording to Eq. [HI and leave the others fixed to their SM 
values. Finally, we choose the (i.e., reconstructed di- 
jet pt) distribution to be our kinematic variable to probe 
anomalous couplings in data. Results are interpreted in 
two different scenarios: LEP parametrization and equal 
couplings, both with A^p = 2 TeV. 



VII. REWEIGHTING METHOD 

The Pythia [32] LO MC generator with CTEQ6L1 
PDFs was used to simulate a sample of WW and WZ 
events at LO. We use the mc@nlo MC generator [33| 
with CTEQ6M PDFs to correct the event kinematics 
for higher order QCD effects by reweighting the differ- 
ential distributions of px(WV) and AR(W, V) produced 
by Pythia to match those produced via mc@nlo. We 
simulate the LO effects of anomalous couplings on the 
Pt distribution by reweighting the SM predictions for 
WW and WZ production from Pythia to include the 
contribution from the presence of anomalous couplings. 
The anomalous coupling contribution to the normaliza- 
tion and to the shape of p|? distribution relative to the 
SM is predicted by the HZW LO MC generator. 

The reweighting method uses the matrix element val- 
ues given by the generator to predict an event rate in 
the presence of anomalous couplings. More precisely, an 
event rate (R) is assigned representing the ratio of the 
differential cross section with anomalous couplings to the 
SM differential cross section. Because the HZW gener- 
ator does not recalculate matrix element values, we use 
high statistics samples to estimate the weight as a func- 
tion of different anomalous couplings. Thus, we consider 
our approach to be a close approximation of an exact 
reweighting method. 

The basis of the reweighting method is that, in general, 
the equation of the differential cross section, which has 
a quadratic dependence on the anomalous couplings, can 
be written as 



da = const ■ \M\ 2 dX 



= const 



\ jvi \smTmW, 



-dX 



\SM 

const ■ \M\ 2 SM [1 + A(X) Ak + B(X)Ak 2 
C(X)\ + D(X)\ 2 + E(X)Ak\ + ...}dX 
da S M ■ R{X; Ak, A, ...) 



(8) 
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where da is the differential cross section that includes 
the contribution from the anomalous couplings, dasM 
is the SM differential cross section, X is a kinematic 
distribution sensitive to the anomalous couplings and 
A(X), B(X), C(X), D(X), and E(X) are reweighting 
coefficients dependent on X. 

In the LEP parametrization, Eq. [8] is parametrized 
with the three couplings Ak 7 , A, and Agf and nine 
reweighting coefficients, A(X) — I{X). Thus, the weight 
R in the LEP parametrization scenario is defined as 

R (X; Ak, A, Agi) = 1 + A(X)An 

+ B(X)(Ak) 2 + C{X)X + D(X)X 2 

+ E{X)Ag 1 +F{X){Ag l ) 2 + G{X)AnX (9) 

+ H(X)AnA gi +I(X)XAg 1 

with Ak = Ak 7 , A = A 7 = Xz, and Agi = Agf. 

In the equal couplings scenario, Eq. [5] is parametrized 
with the two couplings An and A and five reweighting 
coefficients, A(X) — E(X). In this case the weight is 
defined as 

R (A;Ak,A) = 1 + A(X)An + B(X)An 2 . . 
+ C(X)X + D(X)X 2 + E{X)AkX [W) 

with Ak = Ak 7 = Ak^ and A = A 7 = Xz- 

The kinematic variable X is chosen to be the pt of 
the qq system, which is highly sensitive to anomalous 
couplings, as demonstrated in Fig. |4j Depending on the 
number of reweighting coefficients, a system of the same 
number of equations allows us to calculate their values 
for each event. Applied on the SM distribution of X for 
any combination of anomalous couplings, the distribu- 
tion of X weighted by R corresponds to the kinematic 
distribution in the presence of the given non-SM TGC. 

To calculate reweighting coefficients in the LEP 
parametrization scenario, we generate nine different func- 
tions, Fi (i = 1 — 9), fitting the shape of the p q ^ distribu- 
tions in the presence of anomalous couplings. The values 
of anomalous TGCs are chosen to deviate ± 0.5 relative 
to the SM as shown in Table HT1 We calculate nine weights 
Ri normalizing the functions Fi with the cross sections 
given by the HZW generator. 



TABLE II: The values of Ak 7 , A and Apf used to calcu- 
late the reweighting coefficients A(X) — I(X) in the LEP 
parametrization scenario. 





Fi 


F 2 


F 3 


Fi 


F B 


F 6 


Fr 


F 8 


F 9 


Ak 7 








+0.5 


-0.5 








+0.5 


+0.5 





A 


+0.5 


-0.5 














+0.5 





+0.5 


A 5 f 














+0.5 


-0.5 





+0.5 


+0.5 



To verify the derived reweighting parameters, we cal- 
culated the weight R for different Ak, A, and/or Agf 
values, applied the reweighting coefficients and compared 



reweighted pf? shapes to those predicted by the genera- 
tor. Discrepancies in the shape of less than 5% and 
in normalization of less than 0.1% from those predicted 
by the generator represent reasonable agreement. 

When measuring TGCs in the LEP parametrization, 
we vary two of the three couplings at a time, leaving 
the third coupling fixed to its SM value. This gives the 
three two-parameter combinations (Ak, A), (Ak, Agf), 
and (A, Agf). For the equal couplings scenario there 
is only the (Ak, A) combination. In each case, the two 
couplings being evaluated are each varied between -1 and 
+1 in steps of 0.01. For a given pair of anomalous cou- 
pling values, each event in a reconstructed dijet pr bin is 
weighted by the appropriate weight R and all the weights 
are summed in that bin. The observed limits are deter- 
mined from a fit of background and reweighted signal MC 
distributions for different anomalous couplings contribu- 
tions to the observed data using the dijet px distribution 
of candidate events. 



VIII. SYSTEMATIC UNCERTAINTIES 

We consider two general types of systematic uncertain- 
ties. Uncertainties of the first class (type I) are related 
to the overall normalization and efficiencies of the var- 
ious contributing physical processes. The largest con- 
tributing type I uncertainties are those related to the 
accuracy of the theoretical cross section used to normal- 
ize the background processes. These uncertainties are 
considered to arise from Gaussian parent distributions. 
The second class (type II) consists of uncertainties that, 
when propagated through the analysis selection, impact 
the shape of the dijet pr distribution. The dependence of 
the dijet px distribution on these uncertainties is deter- 
mined by varying each parameter by its associated un- 
certainty (±1 s.d.) and reevaluating the shape of the 
dijet pt distribution. The resulting shape dependence is 
considered to arise from a Gaussian parent distribution. 
Although type II uncertainties may also impact efficien- 
cies or normalization, any uncertainty shown to impact 
the shape of the dijet pr distribution is treated as type 
II. Both types of systematic uncertainty are assumed to 
be 100% correlated amongst backgrounds and signals. 
All sources of systematic uncertainty are assumed to be 
mutually independent, and no intercorrelation is propa- 
gated. A list of the systematic uncertainties used in this 
analysis can be found in Table (TTTJ 

IX. ANOMALOUS COUPLING LIMITS 

The fit utilizes the Minuit [34[ software package to min- 
imize a Poisson y 2 with respect to variations to the sys- 
tematic uncertainties [29J]. The \ 2 function used is 




n Ab C p {di]mj{R)) -p- £ G (R k a k ;0, a k ) \ 
^ CPfadk) | = \ r G (0;0,<7 fc ) J 
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TABLE III: Systematic uncertainties in percent for Monte Carlo simulations and multijet estimates. Uncertainties are identical 
for both lepton channels except where otherwise indicated. The nature of the uncertainty, i.e., whether it refers to a normal- 
ization uncertainty (type I) or a shape dependence (type II), is also provided. The values for uncertainties with a shape 
dependence correspond to the maximum amplitude of shape fluctuations in the dijet pt distribution (0 GeV < pr < 300 GeV) 
after ±1 s.d. parameter changes. However, the full shape dependence is included in the calculations. 



Source of systematic 
uncertainty 


Diboson signal 

{%} 


W+jets 
[%] 


Z±jets 
{%] 


Top 

[%} 


Multijet 

[%] 


Tvne 
rype 


Trigger efficiency, electron channel" 


+2/ -3 


±2/ -3 


±2/ -3 


+2/ -3 




I 


Trigger efficiency, muon channel 


+0/-5 


±0/- 5 


+0/-5 


+0/-5 




II 


Lepton identification" 


±4 


±4 


±4 


±4 




I 


Jet identification 


±1 


±1 


±1 


± <1 




II 


Jet energy scale 


±4 


±7 


±5 


±5 




II 


Jet energy resolution 


±3 


±4 


±4 


±4 




I 


Luminosity 


±6.1 


±6.1 


±6.1 


±6.1 




I 


flTOQci Qp.pt~ i ptn finplnH in or PI) n 1 1 ~n pptH~ amtipci 1 




±20 


±6 


±10 




T 


Multijet normalization, electron channel 










±20 


I 


Multijet normalization, muon channel 










±30 


I 


Multijet shape, electron channel 










±7 


II 


Multijet shape, muon channel 










±10 


II 


Diboson signal NLO/LO shape 


±10 










II 


Diboson signal reweighting shape 


±5 










II 


Parton distribution function (acceptance only) 


±1 


±3 


±2 


±2 




II 


ALPGEN r\ and AR corrections 




±1 


±1 






II 


Renormalization and factorization scale 




±1 


±1 






II 


ALPGEN parton-jet matching parameters 




±1 


±1 






II 



a Lepton efficiencies depend on kinematics; however, their frac- 
tional uncertainties are much less kinematically dependent and 
have a negligible effect on the shape of the dijet pj< distribution. 



i=l \ 1 J k=l 

in which the indices i and k run over the number of 
histogram bins (iVj,) and the number of systematic un- 
certainties (N s ), respectively. In this function £ p (a;/3) 
is the Poisson probability for a events with a mean of 
(5 events; C G (x; /i,cr) is the Gaussian probability for x 
events in a distribution with a mean value of /i and a 
variance a 2 ; Rk is a dimensionless parameter describing 
departures in nuisance parameters in units of the associ- 
ated systematic uncertainty er^ ; di is the number of data 
events in bin j- and rrii(R) is the number of predicted 
events in bin i |29j |. 

Systematics are treated as Gaussian-distributed un- 
certainties on the expected numbers of signal and back- 
ground events. The individual background contributions 
are fitted to the data by minimizing this y 2 function over 
the individual systematic uncertainties [29| . The fit com- 
putes the optimal central values for the systematic uncer- 
tainties, while accounting for departures from the nom- 
inal predictions by including a term in the \ 2 function 
that sums the squared deviation of each systematic in 
units normalized by its ±1 s.d. uncertainties. 

Figure O shows the dijet px distributions in the com- 
bined electron and muon channels after the fit. The value 
of x 2 is measured between data and MC dijet px dis- 
tributions as the signal MC is varied in the presence of 
anomalous couplings. The A% 2 values of 1 and 3.84 from 



the minimum x 2 m the parameter space, for which all 
other anomalous couplings are zero, represent the 68% 
confidence level (C.L.) and 95% C.L. limits, respectively. 
For the LEP parametrization, the most probable cou- 
pling values as measured in data with associated uncer- 
tainties at 68% C.L. are k 7 = l-07±g;l|, A = 0.00±° ° 6 6 , 
and gf = 1.041°'^. For the equal couplings scenario the 
most probable coupling values as measured in data with 
associated uncertainties at 68% C.L. are n — 1.041° '\\ 
and A = 0.00±° The observed 95% C.L. limits esti- 
mated from the single parameter fit are -0.44 < Ak 7 < 
0.55, -0.10 < A < 0.11, and -0.12 < A 5 f < 0.20 for the 
LEP parametrization or -0.16 < Ak < 0.23 and -0.11 
< A < 0.11 for the equal couplings scenario (Table ITV|) . 

The observed 68% C.L. and 95% C.L. limits in two- 
parameter space are shown in Figs. [HI and [7] as a function 
of anomalous couplings along with the most probable val- 
ues of Ak, A, and Agf . 

As shown in Table El the 95% C.L. limits on anoma- 
lous couplings Ak 7 , AA, and Agf set using the dijet 
Pt distribution of WW/WZ — > ivjj events are com- 
parable with the 95% C.L. limits set by the DO Col- 
laboration from WW [|, WZ [|, and W-f [5] produc- 
tion in fully leptonic channels using ss 1 fb _1 of data. 
The most recent 95% C.L. one-parameter limits from 
the CDF Collaboration under the equal couplings sce- 
nario at A NP = 1.5 TeV are -0.46 < Ak < 0.39 and 
—0.18 < A < 0.17 using 350 pb _1 of data, combining 
the tvjj and £vf (I = e,fi) final states @. These re- 
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TABLE IV: The most probable values with total uncertainties (statistical and systematic) at 68% C.L. for k 7 , A, and gf along 
with observed 95% C.L. one-parameter limits on A/t 7 , A, and Agrf measured in 1.1 fb -1 of WW/WZ — > Ivjj events with 
A NP = 2 TeV. 



68% C.L. 




A — A-y — Xz 


9? 


LEP parametrization 
Equal couplings 


„ _ 1 n7 +o.26 

l%j — ±.UI _q 29 

Ky = k z = 1.04±g;£ 


A = Q.00±g;gg 
A = Q.0O±g : gg 


9? = 1.04±g;g| 


95% C.L. 


A/t 7 


A A-y — Xz 


Agf 


LEP parametrization 
Equal couplings 


-0.44 < An-, < 0.55 
-0.16 < Ak< 0.23 


-0.10 < A < 0.11 
-0.11 < A < 0.11 


-0.12 < Agf < 0.20 



TABLE V: Comparison of 95% C.L. one-parameter TGC limits between the different channels studied at DO with « 1 fb 1 of 
data: WW -> Ivlv, W7 -» £1/7, WZ -> and WW + -> (i = fM, e) at Ajv P = 2 TeV. 



LEP parametrization 


Ak 7 




A A-y Xz 


Afff 


ITZ-»M (1 fb _1 ) 
Wj^lv-y (0.7 fb _1 ) 
WW -> Ivlv (1 fb _1 ) 
WW + WZ^ ivjj (1.1 fb" 1 ) 


-0.51 < Ak 7 
-0.54 < Ak 7 
-0.44 < Ak 7 


< 0.51 

< 0.83 

< 0.55 


-0.17 < A < 0.21 
-0.12 < A < 0.13 
-0.14 < A < 0.18 
-0.10 < A < 0.11 


-0.14 < Ag? < 0.34 

-0.14 < Ag? < 0.30 
-0.12 < Afff < 0.20 


equal couplings 


Ak 7 




A A-y Xz 


Afff 


WZ^lvll (1 fb" 1 ) 
Wj->£wy (0.7 fb" 1 ) 
WW -> Ivlv (1 ftr 1 ) 
WW + WZ -> (1.1 fb" 1 ) 


-0.12 < Ak 
-0.16 < Ak 


< 0.35 

< 0.23 


-0.17 < A < 0.21 
-0.12 < A < 0.13 
-0.14 < A < 0.18 
-0.11 < A < 0.11 





suits are limited by statistics, but a factor of nearly 10 
times more data is expected to be available for analysis 
by DO by the end of Run II of the Fcrmilab Tevatron. 
With additional data the potential to reach the individ- 
ual LEP2 anomalous TGC limits @, 1, i] shown in Ta- 
ble IVII is significant. The combined LEP2 results still 
represent the world's tightest limits on charged anoma- 
lous couplings [l(| and give the most probable values of 
k 7 , A, and gf as k 7 = 0.973±g;gf| 5 A = -0.028±g;g|5, 
and gf = 0.984±°;°^ a t 68% C.L. 

In summary, we have presented a measurement of 
jWW /ZWW couplings using a sample of semilep- 
tonic decays of WW/WZ boson pairs corresponding to 
1.1 fb -1 of pp collisions collected with the DO detector at 
the Fermilab Tevatron Collider. The measurement is in 
agreement with the SM. On the other hand, this analysis 
yields the most stringent limits on jWW /ZWW anoma- 
lous couplings from the Tevatron to date, complementing 
similar measurements performed in fully leptonic decay 
modes from Wj, WW, and WZ production. 
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TABLE VI: Measured values of k 7 , A and gf couplings and their associated uncertainties at 68% C.L. obtained from the 
one-parameter fits combining data from different topologies and energies at LEP2 experiments. The last column shows the 
DO result obtained from the iv jj final states only selected from 1 fb~ of data. The uncertainties include both statistical and 
systematic sources. 
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FIG. 3: Semileptonic production cross sections for (a) WW and (b) WZ normalized to the SM prediction as a function of 
anomalous coupling Ak (A = Agf — 0) in the LEP parametrization scenario. The new physics scale Ajvp is set to 2 TeV. 




FIG. 4: Normalized distributions of the hadronic W boson (a) pr and (b) rapidity at the parton level in WW production 
including anomalous couplings under the LEP parametrization scenario: Ak 7 = +0.5 (A = Agf = 0, Ak z = -0.15), A = +0.5 
(Ak 7 = Akz = Agf — 0), and Agf = +0.5 (Ak 7 = A = 0, Akz = 1.5) compared to the SM distribution for WW production 
with unity normalization. The new physics scale Ajvp is set to 2 TeV. 
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FIG. 5: (a) The dijet pr distribution of combined (electron+muon) channels for data and SM predictions following the fit of 
MC to data, (b) The difference between data and simulation divided by the uncertainty (statistical and systematic) for the 
dijet pt distribution. Also shown are the MC signals for anomalous couplings corresponding to the 95% C.L. limits for Aft 
and A in the LEP parametrization scenario. The full error bars on the data points reflect the total (statistical and systematic) 
uncertainty, with the ticks indicating the contribution due only to the statistical uncertainty. 
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FIG. 6: The 68% C.L. and 95% C.L. two-parameter limits on the -yWW/ZWW coupling parameters Ak 7 , A, and Agf , in the 
LEP parametrization scenario and Anp = 2 TeV. The dots indicate the most probable values of anomalous couplings from the 
two-parameter combined (electron+muon) fit and the star markers denote the SM prediction. 



13 




Ak 



FIG. 7: The 68% C.L. and 95% C.L. two-parameter limits on the 'yWW/ZWW coupling parameters An and A, in the 
equal couplings scenario and Ajvp —2 TeV. The dot indicates the most probable values of anomalous couplings from the 
two-parameter combined (electron+muon) fit and the star marker denotes the SM prediction. 
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